

LOAM COPY: RETURN TO 
APffL (WLOL) 
KIRTLAND APB, N MEX 


A STUDY OF RADIATION HEAT TRANSFER 
FROM A CYLINDRICAL FIN WITH 
BASE SURFACE INTERACTION 

by Charles A. Cothran 

George C. Marshall Space Flight Center 
Marshall, Ala. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. • APRIL 1970 


0ESTETO 


TECH LIBRARY KAFB, NM 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA TN D-5722 

3. RECIPIENT'S CATALOG NO. 

4. TITLE AND SUBTITLE 

A Study of Radiation Heat Transfer From a Cylindrical 
Fin With Base Surface Interaction 

5. REPORT DATE 

April 1970 

6. PERFORMING ORGANIZATION CODE 

7. AUTHOR (S) 

Charles A. Cothran 

8. PERFORMING ORGANIZATION REPORT # 

M362 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 

Marshall Space Flight Center, Alabama 35812 

Astronautics Laboratory, Science and Engineering Directorate 

10. WORK UNIT NO. 

124-09-05-00-62 

11. CONTRACT OR GRANT NO. 

13. TYPE OF REPORT & PERIOD COVERED 

Technical Note 

12. SPONSORING AGENCY NAME AND AOORESS 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

i ' ' 

14. SPONSORING AGENCY CODE 


Submitted to the University of Alabama in partial fulfillment of the requirements for the 
Master of Science Degree in Engineering, August, 1969 


16. ABSTRACT 

This paper presents a study of the temperature distribution of a cylindrical fin 
transferring heat simultaneously by internal conduction and by surface radiation. The 
temperature distribution was studied considering the fin only and the more complex 
case of surface interactions with a circular base. 

Equations are derived for the temperature distribution of the fin only, the fin 
with base surface interactions when both surfaces are black, and also the case of gray 
diffuse surfaces. Solutions to the equations are obtained by numerical techniques, the 
analysis is verified experimentally, and the results are given. 


17.’ KEY WORDS 

Radiation Heat Transfer 
Base Surface Interaction 
Cylindrical Fin 


18. DISTRIBUTION STATEMENT 

Unclassified - Unlimited 


19. SECURITY CLASSIF. (of thi« r*port> 

20. SECURITY CLASSIF. (of thi. page) 

21. NO. OF PAGES 

22. PRICE * 

Unclassified 

Unclassified 

93 

$ 3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 

Springfield, Virginia 22151 


I 




TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 2 

ANALYSIS 6 

Formulation of Equations for Black Surfaces 6 

Numerical Solution for Black Surfaces 13 

Formulation of Equations for Nonblack Surfaces 13 

Numerical Solution for Nonblack Surfaces 16 

EXPERIMENTAL INVESTIGATION 20 

Base 21 

Fin 22 

Test Setup 23 

Experimental Procedure . 24 

RESULTS 28 

Test Cases 28 

Parameter Study 33 

CONCLUSIONS 35 

APPENDIX A. CONFIGURATION FACTOR ANALYSIS 37 

APPENDIX B. EXPERIMENTAL EQUIPMENT 49 

APPENDIX C. ILLUSTRATIONS 55 


LIST OF REFERENCES 


79 



LIST OF ILLUSTRATIONS 


Figure 

1 . 

2 . 

A-i. 

A-2. 

C-i. 

C-2. 

C-3. 

C-4. 

C-5. 

C-6. 

C-7. 

C-8. 

C-9. 


Title 

Fin System Geometries 

Sketch of Rod-Fin with Circular Base 

Base Surface Area Viewed by Differential Element 
on Rod 

Coordinate Angle Through Which Either dA R Can 

View 6A or dA Can View 6A 

B B -tt 

Test Specimen Mounted in Vacuum Chamber 

Test Specimen 

Schematic of Test Setup 

Tests 1-4 and Calculated Temperature Distributions 
for Rod with Base 

Tests 5-7 and Calculated Temperature Distributions 
for Rod with Base 

Tests 8-13 and Calculated Temperature Distributions 
for Rod with No Base 

Tests 14-17 and Calculated Temperature Distributions 
for Rod with Base 

Temperature Distributions for Thermocouple 

Correction Test Runs, T ( # K) 

Temperature Distributions for Thermocouple 

Correction Test Runs, 8 = T/T 0 


Page 

4 

7 

38 

42 

56 

57 

58 

59 

60 
61 
62 

63 

64 


iv 



LIST OF ILLUSTRATIONS (Concluded) 

Figure Title Page 

C-10. Configuration Factor of a Differential Element 

of Rod Viewing Base from Distance z 65 

C-ll. Temperature Distributions for Fixed N and Various 

Rj and R 2 Ratios 66 

C-12. Temperature Distribution for Fixed N and Various 

Rj and R 2 Ratios 67 

C-13. Temperature Distributions for Fixed N and Various 

Rj and R 2 Ratios 68 

C-14. Some Limiting Cases of Temperature Distribution 

for Black Surfaces ‘ 69 

C-15. Fin Efficiency Variation with Conduction Parameter 

for Fixed Values of Ri and Various R 2 Ratios, 

Rj = 0. 026 70 

C-16. Fin Efficiency Variation with Conduction Parameter 

for Fixed Values of R* and Various R 2 Ratios, 

Rl = 0. 05 71 

C-17. Fin Efficiency Variation with Conduction Parameter 

for Fixed Values of Rj and Various R 2 Ratios, 

Rj = 0. 1 72 

C-18. Fin Efficiency Variation with Conduction Parameter 

for Fixed Values of and Various R 2 Ratios, 

R = 0.5 . . 73 

C-19. Fin Efficiency Variation with Conduction Parameter 

for Fixed Values of R t and Various R 2 Ratios, 

R 1 = 1.0 74 


v 



LIST OF TABLES 


Table Title Page 

C-i. Calculated Temperature Distributions 75 

C-2. Test Summary 76 

C-3. Test Data 77 

C-4. Temperature Error at the Insulated End of the 

Rod 78 


vi 



NOMENCLATURE 


A 

A* 

6A 
B 
B* 
d( ) 

F 
F* 
g( ) 

H 

H* 

A A A 

i, j ,k 
k 

1, m, n 
L 
M 
N 

c 

N 

N* 

A 

n 


Area (cm 2 or ft 2 ) 

That portion of a larger area that is viewed from a specified 
location (cm 2 or ft 2 ) 

Area with differential width and finite length 
Radiosity (W/cm 2 or Btu/ft 2 -hr) 

Radiosity (dimensionless) 

Differential 


Configuration factor (dimensionless) 

Configuration factor in terms of dimensionless parameters 
Irradiation from surroundings (W/cm 2 or Btu/ft 2 -hr) 
Irradiation from the system (W/cm 2 or Btu/ft 2 -hr) 
Irradiation from the system (dimensionless) 

Unit vectors 

Thermal conductivity (W/cm-°K or Btu/ft-hr- 0 R) 
Direction cosines 
Length (cm) 

Number of isothermal segments of the specimen 
Conduction parameter (dimensionless) 



Unit vector normal to a surface 


vii 



I 


q 

R 

r 

S 

A 

S 

s 

T 

X 


x,y,z 


z 

a 

6 .. 

ij 

e 

V 

e 

7T 

P 

P* 

a 


T> T 


NOMENCLATURE (Continued) 

RatB of heat flow (W) 

Dimensionless ratio, r/L 
Radius (cm) 

Line 

Vector 

Magnitude of line S 

Absolute temperature (degrees Kelvin, e K) 

6.. - (1 - e.) 

Matrix of — — F XA XA for 1 < i < M 

e. 6A. - 6A. 

i i 3 

Rectangular coordinates 
Dimensionless ratio, z/L 
Absorptivity (dimensionless) 

Kronecker delta 
Emissivity (dimensionless) 

Efficiency (dimensionless) 

Dimensionless temperature, T/T 0 
3. 14159 

Cylindrical radius coordinate 
Reflectivity (dimensionless) 

Stefan- Boltzmann constant (5.6699 x 10“ 12 W/cm 2 -‘K 4 or 
0. 1713 x 10- 8 Btu/ ft 2 -hr- 6 R 4 ) 

Dummy integration variables 


viii 



I 


NOMENCLATURE (Concluded) 


<p 

4 > 

Subscripts: 

B 

i,j 

R 

<5A 

dA 

0, i, 2. . il 


Cylindrical coordinate angle 
Inverse of X 

Base surface 

i denotes a specific element of j elements 
Rod surface 
Differential strip area 
Differential area 

Thermocouple station locations on the rod test specimen 


ix 


L 


ACKNOWLEDGMENT 

The author would like to take this opportunity to acknowledge the 
contributions of those persons whose assistance was most prominent in the 
preparation of this thesis. Sincere thanks are extended to Dr. A. R. Shouman 
for his many helpful comments and for his efforts in sponsoring the use of the 
MSFC Propulsion and Vehicle Engineering Laboratory test equipment and 
facilities. Sincere thanks are also extended to Dr. R. Blanton, his thesis 
advisor, for his many helpful comments and his guidance and to Miss Lucy 
Coward of the MSFC Computation Laboratory who programmed the numerical 
solutions of the temperature distribution equations. 


x 



A STUDY OF RADIATION HEAT TRANSFER FROM A 
CYLINDRICAL FIN WITH BASE SURFACE 
INTERACTION 

SUMMARY 

In this study the temperature distribution and heat transfer rates of a 
cylindrical radiating fin were investigated, considering no base surface inter- 
action and also the case of interaction with a circular base. Surroundings of 
the fin were assumed to be black and at a temperature of absolute zero. 

The equations governing the temperature distribution of a fin attached 
to a circular base causing radiation interaction were derived for two cases: 

f ( 

when their surfaces were black; and when their surfaces were gray and diffuse. 
Solutions to the equations were obtained by numerical techniques and the 
analysis was qualitatively verified with experimental studies. 

Numerical solutions were obtained for the specific conditions of each 
of the 20 test runs and for a wide variation of the governing parameters for a 
black rod and base. The parameter study showed that the rod heat transfer 
efficiency may vary as much as 50% with changing base size. The experimental 
data and numerical solutions showed trend agreement. However, it was shown 
that substantial error in the measured temperature distribution could be 
attributed to the use of an excessive number of thermocouples and the carbon 
coating used during ten of the tests. 



INTRODUCTION 


The subject of heat transfer from fins and extended surfaces has been 
studied analytically and experimentally for almost two centuries [1] . Until 
recent years most work was concentrated on convection from fins with constant 
thermal properties. The problem of a fin exchanging heat with its surround- 
ings by radiation has only recently come under extensive study as a result of 
man's interest in space and the newly found capability for space travel. 

The temperature distribution along a thin rod, wire, or thin-walled 
tube heated by an electric current while in a vacuum has been of interest, 
both experimentally and theoretically, for a long time. The literature on this 
problem contains work by Langmuir [2], Langmuir and Taylor [3] , Worthing 
[4], Worthing and Holliday [5], Stead [6] , Bush and Gould [7], Baerwald [8], 
Jain and Krishman [9], and Shouman [10]. Shouman [10] presented an exact 
solution for a one-dimensional fin with constant cross sectional area and 
constant thermal properties that transferred heat from its surface by radiation 
only to an equivalent sink temperature. 

The interest in space exploration generated an interest in fins of 
various geometries applicable to heat exchanger designs. The early 
published work in this area treated the temperature distribution of isolated 
fins. Lieblien [11] analyzed the temperature distribution along a rectangular 
fin and Chambers and Somers [12] presented a similar analysis of a circular 
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fin. Tatom [13] presented the mathematical formulation for the temperature 
distribution along a rod and presented an exact solution for the special case 
of an infinitely long rod. 

The analysis of fins then proceeded to the consideration of mutual 
irradiation of fins and their base surfaces. ' Callinan and Berggren [14] con- 
sidered a wide range of design problems and studied the fin geometry shown 
in Figure 1A; however, the base surface effect and multiple reflections were 
not considered. Bartas and Sellers [15] indicated the importance of radiation 
interaction for this design and considered a singly reflected ray in the analytical 
formulation but omitted the singly reflected ray for the solutions presented. 

The problem of flat rectangular fins sharing a common edge, Figure 
IB, was studied by Sparrow, Eckert, and Irvine [16] . These authors formu- 
lated the solutions for this geometry, and in another study [17] , they presented 
solutions for mutual irradiation between fins whose surfaces are gray and 
reflect diffusely. A similar arrangement, but of trapezoidal cross section, 
was considered by Karlekar and Chao [18]; they presented an optimization 
procedure for achieving maximum dissipation. Hering [19] studied this fin 
system and presented solutions for gray surfaces which reflect specularly, 
and he compared the effectiveness of the cases of specular and diffuse 
reflections. 
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(C) 


FIGURE 1. FIN SYSTEM GEOMETRIES 
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A complete mathematical formulation of the governing equations for the 

system shown in Figure 1A was presented by Sparrow and Eckert [20] . They 

considered the mutual irradiation between the fin and base surfaces and 

irradiation from external sources. The case of selective gray surfaces was 

also considered. However, the solutions presented were for black surfaces 

with negligible irradiation from external sources. The physical parameters 

were varied over a wide range and the conduction parameter, N , was varied 

c 

from 0. 5 to 20. Sarabia and Hitchcock [21] presented the solutions for this 

system in terms of overall fin-tube efficiency for the black and gray diffuse 

cases. The conduction parameter was varied from i to 20 for the black case 

and the data were presented for N =1 for the gray diffuse case. Hwang- Bo 

c 

[22] also considered this sytem and presented a method of maximizing heat 
rejection per unit weight; however, irradiation was simplified by neglecting 
components reflected twice or more. 

Sparrow, Miller and Jonsson [23] studied the fin system shown in 
Figure 1C with gray diffuse surfaces. The system consisted of annular fins 
mounted on a constant temperature tube with mutual irradiation between the 
fins and the tube. The problem was solved for a wide range of physical 
parameters, and the conduction parameter was varied from 0 to 0. 7. 

The effects of convection were considered negligible in the previously 
mentioned papers and were neglected for this study. Gordon [24] substantiated 
that convection was negligible at pressures below 10“ 4 millimeters of mercury. 
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The mathematical formulation of the energy balance on a differential 
element of a fin transferring heat by internal conduction and by surface 
radiation results in a second order nonlinear differential equation. The 
equation is further complicated when the mutual irradiation by the fin base 
surface, other fins, and external sources are considered. As a result of these 
complications, an exact solution is not known and numerical techniques must 
be employed. An exact solution for a fin not considering interaction with other 
surfaces has been published by Shouman [10] . In a test program to verify 
his solution for the temperature distribution of a fin, Shouman used an 
aluminum rod 2. 54 cm (1 inch) in diameter and 48. 90 cm (19^ inches) long. 

As an extension of this work it was decided to investigate, both analytically and 
experimentally, the effect of adding a base surface upon the temperature 
distribution of the radiating rod. 


ANALYSIS 

Formulation of Equations for Black Surfaces 

The following analysis considered the heat transferred from a black 
circular rod including interaction with a black base surface. The base surface 
considered was circular with a uniform temperature. The temperature 
gradient in the fin was assumed extremely small in the radial direction com- 
pared to the gradient in the z direction which reduced the problem to one- 
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dimensional heat flow in the rod. In addition, the following assumptions were 


made: 

1. The rod was isotropic and homogeneous. 

2. The material thermal properties were constant. 

3. There were no heat sources or sinks in the rod. 

4. The system and surroundings were at thermal equilibrium. 

Now consider a differential slice of the rod shown in Figure 2A and C 
at a distance z from the base. The net heat, q, conducted into the element dz 




(C) 

FIGURE 2. SKETCH OF ROD-FIN WITH CIRCULAR BASE 
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was equated to the heat lost by radiation, or 


( dq conduction) XT . ( dq radiation) 

Net Net 


(II-l) 


The amount of heat that entered the element at z by conduction was 


- n aT 

z) = - r R ( k T^ 


( H-2) 


where T represented the absolute temperature and k, the thermal conductivity. 

To determine the corresponding quantity of heat conducted from the element 

9T 


at z + dz, let q( z, T) = k • 


9z 


and allow z to increase by dz. Then when 


q(z + dz, T) was expanded in a Taylor's series and all terms after the first 
two were neglected [25] , 


q(z + dz, T) = q(z, T) + dz 


(H-3) 


or 


. di, T| =.-£-+ ^ 


(•f-) 


dz 


(11-4) 


Then, the heat conducted at z + dz was 


, 2 r, 9T 3 /. 9T \ , 

da = - 7r r^ k 1- — ( k ) dz 

q (z + dz) R [ 9z 9z V, 3z / 


( II— 5) 
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and the net conduction into the element was 


(^conduction) Net “ ' r R ( k l£) d2 ] ( “- 6) 

or 

2 

(dq , .. ) =7rr2k^dz (II-7) 

V conduction y . R 3z^ 

since k was assumed constant. 

When the radiation from the fin was considered, there were three 
possibilities: 

1. The surface of the element emitted a quantity of radiant energy 

ex T 4 (27rr D dz), (II-8) 

where cr was the Stefan- Bo ltzm an constant. 

2. There was a quantity of energy incident on the element due to the 
base surface. When the configuration factor, F, which denotes 
the fraction of energy leaving one surface which was incident on 


another was introduced, the base surface contribution became 



and when the reciprocity relationship between configuration 
factors, 


A_F. . =27rr„dzF ! .. 

B a b - 6A r r 6A r 


l B’ 


(II- 10) 


was used the energy incident on the element of rod from the base 
was written as 


-^ E (2^ R dz, F • (II-l 1) 

R B 

3. There was a quantity of energy from some external source. 

Since the base surface would probably block a portion of any energy 
from an external source, the incident energy would vary with z. 
When the external radiation was defined as g(z) per unit area, 
the quantity incident on the element dz was 


-g(z) 27rr R dz. 


( 11 - 12 ) 


The net radiation transfer was written as 


(^radiation) 


Net 


R B 


27rr dz. (11-13) 
R 


When combined, the net conduction and the net radiation transfer yielded 


JEv J3 


(H-14) 


which was written as 
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(11-15) 


d 2 T _2o_ ^ 1 2g(z) 

^ r R k L B 6A r- a bJ r R k 

The incident energy from the surroundings may be neglected if the rod 
temperature is sufficiently high compared to the surrounding temperature. 
Since the tests were conducted with the surroundings at the temperature of 
liquid nitrogen, the term for incident energy was neglected. 

The differential equation for the temperature distribution in the rod 
was written in dimensionless form by introducing the following nondimensional 
quantities, 


T T B „ z r R , r r B 

’‘T.' 8 b = t 7 ' Z ’T- R >=— and r ’ = T" 


Substitution into equation (11-15) and neglect of g(z) yielded 


d 2 (Tf,e) 2a 
d(ZL) 2 r R 


k [(T,9)*-(T 0 9 b )*F* A r _ a J 


(11-16) 


which reduced equation (11-16) to 


d z e _ 2gL 2 T n 3 


dZ 


'fll r fl 4 

k L " 


R 
2rp 3 


B 6A r -A b 


] 


(H-17) 


2crLrT 0 

The quantity r-°- was defined as N and is generally referred to as the 

r R k ° 

conduction parameter in the literature. The above equation then became 


d 2 0 

dZ? 


= N o [«,* - 


fl 4 IT* 

B 6A R- A B 


] 


(11-18) 


The term 6 was equal to 1 at Z = 0 when the rod temperature was equal to 
B 

the base temperature. For the tests conducted, however, 0 was not equal 

B 


li 



to 1. The derivation of the configuration factor F* was quite lengthy 

6A r “ a b 

and is presented in detail in Appendix A. 

The expression for the temperature distribution from a black rod 

viewing a black circular base was completed by defining the boundary con- 

dT 

ditions: at z = 0 the temperature T = T 0 and at z = L, -j^- = 0 and 

dfl 

correspondingly 0 = 1 at Z = 0 and = 0 at Z = 1. 

The heat transferred from the rod was found quite simply since at 
thermal equilibrium the heat loss from the rod surface must be equal to the 
amount of heat conducted into the rod at the base, or 


Conduction 


(^radiation) 


= q 


Net 


R 


or 


( 11 - 19 ) 


q E =-kU^) 



( 11 - 20 ) 


When the above equation was written in terms of the dimensionless variables 
6, Z, and N , the expression for q R became 

q R i /d6 \ 

oT{(2 T r R L) '-ir(dz ) z=0 (U - 

Examination of the left-hand side of the above equation revealed that the 


quantity was the ratio of the actual fin heat transferred to the heat transferred 
from an ideal fin which was the fin efficiency, 17. 



Numerical Solution for Black Surfaces 


The equation (II-18) for the temperature distribution was solved 

numerically using the Runge-Kutta method. It was solved by determining 

d0 

which of a succession of guessed values of — — at Z = 0 led to a solution 

Q Zu 

d0 

that satisfied the boundary condition = 0, at Z = 1. A AZ step size 
of 0. 01 was used for all solutions; however, AZ step sizes of 0. 001 and 
0. 0001 were tried and it was determined that die accuracy was not improved. 

Solutions were determined for each of the test cases and also for a wide selec- 

r r 

T? R 

tion of N and the dimensionless ratios, R< = — — and R? = — . which 

c 1 L L 

governed the configuration factor, F A * ^^ ese data are presented in 

R ~ A B 

the curves in Appendix C and typical values of Q( Z) are presented in Table 
C-l of Appendix C. 

Since the Runge-Kutta integration of equation (11-18) yielded , 

z 7 Z=0 

the fin efficiency (77) or the fin heat transfer (q ) was easily obtained from 
equation (11-21) . The fin efficiency was plotted versus the conduction param- 
eter for the limiting cases of R 2 = R A (no base) , R 2 = 00 (infinite base, 

F . =0.5) and two intermediate values of R 2 for each of five values 

6A r ~ a b 

of R t . The results are plotted in Appendix C. 


Formulation of Equations for Nonblack Surface 


In this section the analysis is extended from the black case to the case 

of gray, diffuse emitters and absorbers. When a surface element on the rod 

a distance z from the base is considered, the energy that leaves this location 
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includes both the direct emission and the reflected incident energy. This is 
called the radiosity B (energy/ unit area-time) , 


B = eo-Ti + p'H 

ri 


(11-22) 


where H is the incident energy per unit time, and e and p f are the 
emissivity and reflectivity, respectively. 

The net heat loss from a surface location was B - H and the fin energy 
balance was, 


d 2 T 

d?" 


= IT r (b h 


V 


(11-23) 


where the conduction was evaluated the same as in the black analysis. T and 
Z were normalized as before 


6 = and Z = 


which yielded 


d*6 


2 1 / 


dZ* kr R T 0 


< b r-V 


(11-24) 


When equation (11-24) was multiplied and divided by aT^, and B 1 ' and 1 1 ' 

B H 

were defined as B* = — rj- and H* = — — r- , then 

cr 1 o o i 0 


d 2 g 

dZ 


^ =N c (B i- H R } 


(11-25) 


The irradiation H* was a function of the base surface radiosity and of the 
li 

radiation from the surroundings; thus 
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(11-26) 


+ i b z F * 


6A^ - 6A 
-K 


B 


When selective surfaces were not considered, then 




(11-27) 


where e and p were gray body properties, that is, 


e = a = i - p' 


(11-28) 


The radiosity distribution B* in equation (11-26) was related to the other 

13 

energy transfers of the system by the following radiant flux balance 


b b- £ b 9 b + p ' 


[ B* F* 
O-TJ R 6, 


R 


R M b -6A r 


(11-29) 


The four equations (11-25), (11-26), (11-27), and (11-29) in four unknowns 

described the temperature distribution of the rod. However, H* was 

R 

eliminated by use of equation (H-27). 


R P R 


(11-30) 


Substitution into equation (11-25) yielded 

°R 


d 2 e 

dZ 




( II— 3 i ) 


and similarly equation (11-26) yielded 


P'j, g( z ) , 

B H = e R 04 + oTg + P R I B B F 6A - 6 A 

A B R 1 


(11-32) 
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This yielded three equations in three unknowns, B* B* , and 6 ( Z) , the 

-tv B 

temperature distribution, which were to be determined. 

The problem was completed by defining the boundary conditions. At 

dT 

z = 0 the temperature T = T 0 , and at z = L, — — = 0. Correspondingly, 

QZ 

6 = 1 at Z = 0, and — zr = 0 at Z = 1. This nonlinear system of differential 

Q La 

and integral equations required a simultaneous solution. 

Since the case being considered had negligible incident radiant energy, 
g(Z) « 0, the system of equations was reduced to 



(11-33) 



The above system of equations could have been solved numerically in 

this form by assuming distributions for B* and B* and iterating until the 

rv B 

boundary conditions were satisfied. This would have been quite tedious, 
however, since equations (11-34) and (11-35) must first be satisfied. Then, 
BZ' must be substituted into equation (11-33) to find 0( Z) which must be used 


as an input into equations (11-34) and (11-35) until agreement is reached. 
16 



This process must be repeated until the boundary conditions are satisfied. 
Because of the lengthy iteration involved in solving the system in this form, 
it was more expeditious to convert equation (11-33) into an integral form 
which included the boundary conditions. 

When N was defined as 


N = 


N e T5 

c R 


P'x 


(11-36) 


R 


equation (11-33) became 


n[«‘(Z) -B*(Z)] . 


(11-37) 


The above equation was integrated from 1 to Z and the boundary condition 

) = 0 applied to yield 

\dZ/ „ . 


Z=1 


dfl(Z) 

dZ 


z, 

= N f |V(r) - B*(r)] dr 


(11-38) 


where r was a dummy integration variable. Equation (11-38) was integrated 
from 0 to Z and the boundary condition 6=1 at Z = 0 was applied to yield 

Z r' 


0( Z) = 1 + N f f [flV) - B*(t) | dr dr'. 


t (T) . 


(11-39) 


The above double integral was transformed into two first order integrals using 
integration by parts as outlined by Hering [19] to yield 

Z 

0 


6(Z) = 1 + N f (Z-r) |V(t) - B*(r)J dr 
- f Z [fl 4 (r) - Bj(r)] dr . 


(11-40) 
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The above equation could be readily solved by numerical techniques if B^(t) 
were known. 

The procedure outlined by Sparrow and Cess [26] was used to 

determine B*(t) . The rod and base were considered to be M isothermal 
R 

elements. When the_ith element was considered, the integral equation (11-35) 

for B*(r) was expressed as a summation 
R 

j=M 

B? = e. 0 4 +p\ V B* F (H-41 

l il l/-' J 6A. - 6A. 

3=1 i J 

where l<isM, Equation (11-41) was rewritten as 
j=M 

y x. b: s = 0 4 (h-42 

3=1 13 1 1 


where 


X.. = 
i] 


6 ij (1 V F 6A. - 6A. 

J i l 


The term 6 is the Kronecker delta which is 1 when i = j, and zero when 
ij 

i # j. The coefficients, X„, of equation (11-42) form an M- by -M matrix 
where 1 < i < M and 1 £ j ^ M 


x 

X 

x 

11 

12 

1M 

X 

X 

... x„„ 

21 

22 

2M 


X M1 X M2 X MM 
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The inverse matrix of X was defined as tp 9 




'li 

*12 ••• 

” *1M 



'21 

^22 * * * 

“ ^2M 


(11-45) 

'mi 

^M2 '*• 

’ * ^MM 

. 



The coefficients, ip^ 9 were used and equation (11-42) was rewritten as 


j=M 

b* = V i p..e i 

i li i 

j=i J J 


(11-46) 


The above equation was used to determine B* ( t) of equation ( 11-40) . The 

Jti 

specimen was broken into 11 isothermal elements along the rod and 9 on the 
base surface. This technique for determining B*(r) allowed equation (11-40) 
to be numerically integrated to yield the temperature distribution along the 
rod. 


No attempt was made herein to elaborate the above technique for 
calculating B * , and only the steps needed to complete the computer program 
were presented. 

Equation (11-40) was numerically integrated using the trapezoidal 
rule where the rod was broken into 11 segments and the base, 9. An initial 
value of 0 = 1 was assumed to calculate B* . Then B* was used to calculate 
a 0 value which was used as a second input. The process was repeated until 
agreement was reached. 
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Since this solution did not directly yield the value , as did 


Z=1 


the Runge-Kutta solution for the black case, the heat transfer was calculated 
by determining the slope at Z=0 using the following numerical differentiation 
formula [27]: 





(11-47) 


The subscripts refer to the beginning points of the first five isothermal 
elements of the rod. By using this value for , the rod heat 

transfer was determined as before, using equation (11-21) . 


EXPERIMENTAL INVESTIGATION 


The experimental setup consisted of the test specimen, the high 
altitude simulation system, and the data acquisition system as illustrated in 
Figure C-3 of Appendix C. A detailed description of the high altitude simula- 
tion system and the data acquisition system is presented in Appendix B. 

The test specimen consisted of two prime parts — the 2. 54-cm (1-inch) 
diameter aluminum rod and the flat circular 20. 32-cm (8-inch) diameter 
aluminum base. Both the base and the rod had individual heating elements. 

The base was designed in two halves so that it could be mounted on the rod 
without disturbing the heater and the thermocouples previously used to 
determine the temperature distribution of the rod alone. 
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Base 


The base surface consisted of two halves of a circular plate (mounting 
attachments so that the base could be fixed to the rod) , a heater element on 
each half, and 12 thermocouples spot-welded to the base surface viewing the 
fin as shown in Figure C-2. Three different heater arrangements were 
bench tested before a satisfactory temperature distribution was obtained. The 
configuration selected was one with heater wires closely spaced in a circular 
pattern. 

The base heater was constructed in the following manner. Two heater 
wire terminal posts were affixed to each half of the base and the surface was 
then covered with a layer of fiberglass tape to electrically insulate the heater 
wire and base surface. Then 16-gauge Chromel-A resistance wire with a 
resistance of 0. 830 ohms per meter (0. 253 ohms per foot) was wound 
circumferentially on each half of the base. The heater wire was covered by 
a layer of asbestos (0. 318 cm thick) ; this was held in place by a piece of 
aluminum sheet stock (0. 079 cm thick) attached to the base by screws. The 
resistances of the base surface heaters were measured and found to be 5 ohms 
each. 

For bench testing the base specimen heaters were connected to a 
30-volt, 10-amp, DC power supply. The current to the heater wire was 
controlled by a variable resistance. The thermocouple cold junction reference 
was an ice bath and the thermocouple millivolt output was recorded on a strip 
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chart by a Brown-Marco recording system. Several tests were conducted at 
different power levels and the base temperature spread was determined after 
steady state was reached. Steady state was normally achieved in about three- 
fourths of an hour. The tests were conducted at room temperature of 
approximately 297. 6°K (76° F) . 

The base heater configuration selected produced a maximum tempera- 
ture differential of 8°F at an average surface temperature of 388. 7° K (240® F). 
This was considered acceptable since the temperature distribution was expected 
to be much more uniform in vacuum. 


Fin 

The rod or fin portion of the test specimen consisted of an aluminum 
alloy rod, a heater element, and 48 thermocouples as shown in Figure C-i. 

The rod specimen was a 2.54-cm (1-inch) diameter 6061 T6 aluminum alloy rod 
approximately 60.96 cm (24 inches) long. One end of the rod was circum- 
ferentially wound for approximately 10. 16 cm (4 inches) with a heater wire; the 
remainder of the rod contained 48 thermocouples evenly spaced at 12 stations; 
the end of the rod opposite the heater had 4 thermocouples spot-welded to the 
end which was covered with one-fourth inch of insulation. The 48 thermocouples 

3 

were located at intervals of 4.45 cm (1— inches) with 4 couples at each station. 
The 4 thermocouples at each station were evenly spaced around the circumference 
of the rod. 

The fin heater was constructed in the following manner. First, the rod 
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was wrapped with a layer of fiberglass tape to electrically insulate the heater 
wire from the rod. The heater wire was then attached to one terminal post 
and wound around the rod and attached to the second terminal post. It was then 
coated with a layer of Ceramacoat 512 (AREMCO Products, Inc. , Briarcliff 
Manor, New York) and wrapped with a 1. 27-cm (1/2-inch) layer of asbestos 
cloth. The heater resistance was measured and found to be 1. 16 ohms. 

The power input to the rod and also to the base was a closely regulated 
DC power supply. Variable resistances were placed in both base and rod 
power leads so that the desired power supply could be accurately achieved. 

The power input was determined by two Hewlett-Packard model 3430A digital 
voltmeters. One measured the system voltage drop and the other measured the 
voltage across a known resistance placed in one leg of the circuit. 

Test Setup 

The test specimen was suspended in the vacuum chamber by a cord at 
either end of the specimen which attached to a horizontal cross member as 
shown in Figure C-i. The test specimen thermocouple leads were taped to the 
supporting cross member and then fed through the chamber wall with a standard 
vacuum feedthrough connector. A thermocouple cold junction, ice bath, was 
located adjacent to the chamber and then the thermocouples were connected with 
the data acquisition system. The test specimen power leads were similarly fed 
through the vacuum chamber wall and connected to the regulated DC power 
source. Separate vacuum feedthrough connectors were used for the power 
leads and thermocouples to avoid any interference. 
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The vacuum chamber walls were painted black to minimize the reflected 
energy. To prevent the complications of having to account for an incident 
energy distribution on the specimen, the vacuum chamber wall was maintained 
at 78. 7°K (-318° F) by circulating liquid nitrogen through the vacuum chamber 
wall. The temperature distribution of the rod without the base was solved 
numerically, including the sink temperature of 78. 7° K to verify that the 
incident energy was negligible. The inclusion of the incident energy was 
calculated to produce less than 0. 1% change in the temperature at the end of 
the rod. 


Experimental Procedure 

The test was designed to explore the temperature range from ambient 
to approximately 533° K (500° F). Based on the eariler rod tests by Shouman, 
approximate values of the rod temperature were known for various power inputs 
to the rod heater. Also, it was known that approximately 8 hours were 
required to achieve steady state heat transfer. Seven of the 20 test runs were 
conducted on a single 8-hour shift. For these runs the rod and base power 
supplies were left on overnight to permit faster achievement of steady state 
conditions. Steady state was considered achieved when all temperatures on 
the rod varied less than 1 degree per hour. This assumption was well borne 
out by longer runs that permitted a steady state condition where essentially 
no temperature variations were noticeable. 
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The test runs are summarized in Tables C-2 and C-3 of Appendix C. 

Test runs 1 through 7 were conducted with the base installed. Test runs 8, 9, 
and 10 were conducted with the rod only, and runs 11, 12, and 13 were con- 
ducted with the rod only and a smoke black carbon coating. Test runs 14 
through 17 were conducted with the base installed and both the base and the 
rod carbon coated. Test runs 18, 19, and 20 were a repeat of run 12 to study 
the effect of removing the thermocouples. All test runs, except the first one, 
were made with liquid nitrogen cooling the chamber walls. 

The initial test was conducted without liquid nitrogen and without 
power to the base to observe the effect of the presence of the base. This was 
an abbreviated test run and the thermal equilibrium was questionable. Test 2 
was a repeat of test 1 but with liquid nitrogen in the cooling jacket. 

Since the test specimen had individual heaters on the rod and base, it 
was highly impractical from a time standpoint to achieve equal temperatures 
on the base and the fin at the base fin intersection. Because of this problem 
the tests with the base were planned to achieve a somewhat higher and a some- 
what lower temperature than the fin temperature at the base. Test 3 and test 4 
satisfied this requirement for the initial power setting to the rod. For tests 5, 
6, and 7 the procedure for tests 2, 3, and 4 was repeated with a different power 
supply to the rod heater as shown in Table C-2 of Appendix C. 

The base was removed and tests 8, 9, and 10 were conducted with 
three different power inputs to the rod heater. These data were desired so 
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that the case of base surface interaction could be compared to the rod only, 
without interactions. 

A carbon black coating was applied to the rod and base to obtain data 
that would satisfy a black surface analysis. The coating was applied by holding 
an acetylene torch burning an acetylene rich mixture near the specimen and 
rotating it until a velvety black coating was obtained over the entire surface of 
the rod and base. Tests 11, 12, and 13 were then made using the rod only. 
These tests were similar to tests 8, 9, and 10. Tests 14, 15, 16, and 17 were 
then conducted with the base installed. The procedure was again to obtain data 
with the base temperature somewhat above and below the rod temperature. 

Test runs 18, 19, and 20 were then conducted with the rod only to 
ascertain the effect of the thermocouples mounted on the rod. The power 
input for these tests was maintained the same as for test 12. For test 18, 
one-half of the thermocouples were removed leaving 2 thermocouples at each 
of 12 stations. For test 19, one-half of the remaining thermocouples were 
removed leaving 1 thermocouple at each station except the fifth station which 
had none. Test 20 was conducted with 1 thermocouple at the heater end and 
1 at the insulated end of the rod. 

The test procedure was essentially the same for all test runs once the 
specimen was prepared for that particular test. The test runs were begun by 
simultaneously turning on the data acquisition system, the power supply to the 
specimen heaters, the backing pump, the heaters to the ejector pump and 
diffusion pump, the cooling water to the pump condensers, and then pressurizing 
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the liquid nitrogen storage trailer (for test runs 1 through 7 the specimen 
heater power supply was turned on the previous day and left on all night) . 

When the diffusion and ejector pump oil reached approximately 545® K (520® F), 
the roughing pump was turned on and the liquid nitrogen flow to the pumping 
system cold trap and to the vacuum chamber cooling jacket was begun. The 
vacuum was then read on the Pirani gauge and the system was switched to 
"hi-vacuum, " or onto the diffusion and ejector pumps when the vacuum reached 
150 microns. 

After hi-vacuum operation was achieved, the procedure consisted of 
checking the vacuum, thermocouple ice bath, heater power supply readings, 
and the specimen temperature distribution periodically. The specimen 
temperature was observed until steady state was achieved and then a new test 
run was begun by changing the power inputs to the specimen. Once a test run 
or series of runs was completed all of the test equipment was turned off and 
the nitrogen storage trailer vented. The cooling water to the ejector and 
diffusion pump oil condenser was left on for approximately 2 hours to protect 
solder joints. 
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RESULTS 


Test Cases 

Twenty tests were conducted as described in the experimental chapter, 
and equations (11-18) and (11-40) were solved (using the rod temperature at the 
base location as the initial temperature value) to determine the correlation 
between the analytical predictions of the temperature distributions and the 
experimental data. The test and calculated temperature distributions were 
plotted in Figures C-4 through C-9 of Appendix C. Since the fin root tem- 
perature and base temperature values for the test cases were used as starting 
points for the analytical solutions, the error between the test data and the 
calculated temperatures was greatest at the end of the fin and was used as a 
criterion to determine agreement between the two. The test data agreement 
with the analytical predictions ranged from -2. 34 % to +5. 97 % as shown in 
Table C-4 of Appendix C. 

By examination of the percentage error for each test as presented in 
Table C-4 of Appendix C, it was seen that the first 10 tests showed a smaller 

error than the last 10 and that the error for the first 10 tests was negative 
and was positive for the last 10. The first 10 tests were conducted with the 

specimen uncoated and the last 10 tests were conducted with the specimen 
coated. Also it was shown by tests 12, 18, 19, and 20 (Figures C-8 and C-9 
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of Appendix C) that the presence of the thermocouples affected the temperature 
distribution. The coating effect and the thermocouple effect were then in- 
vestigated further. 

The thermocouple effect was experimentally investigated in the last 
3 test runs by maintaining the power input to the rod constant at the value used 
in test 12. The thermocouples were progressively removed as described in 
the experimental chapter. Examination of the temperature distribution curves 
of Figures C-8 and C-9 of Appendix C showed that the rod temperature was 
forced higher to dissipate the constant heat input when the thermocouples were 
removed. The fin temperature increased approximately 6. 7° K (12° R) at the 
root of the fin and approximately 3. 9°K (7°R) at the insulated end when the 
48 thermocouples were reduced to only 1 thermocouple at each end of the rod. 

The heat transfer was then calculated from equation (11-21) using the 
slope that was determined by the numerical solution of equation (11-18) . The 
calculated heat transfer showed a progressive increase with the removal of 
the thermocouples even though the heat input to the specimen was constant. 

This was expected since the math model did not compensate for the thermo- 
couple and an increase in the rod base temperature would result in a greater 
calculated initial slope. This indicated increase in heat transfer was 
attributed directly to the thermocouples. The heat transfer from the rod for 
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the thermocouple correction test cases based on the calculated temperature 


distributions was: 


Test 

q R (W) 

q R (Btu/hr) 

12 

19. 77 

67. 5 

18 

21.47 

73. 3 

19 

22. 00 

75. 1 

20 

22. 14 

75. 6 


These results showed that a 10. 7 % change in heat transfer resulted from 
the presence of the thermocouples. 

The test data were then used to determine the heat transfer for the 
thermocouple test runs. The initial slopes for test runs 12, 18, and 19 
were determined by the numerical differentiation formula (11-47) . Even though 
the heat input was constant to the rod, the slope was expected to vary for each 
test since the heat loss at each thermocouple station would decrease as the 
thermocouples were removed. Equation (11-21) was then used to calculate 
the heat transfer from the rod which was: 


Test q R (W) q R (Btu/hr) 


12 

13. 36 

45. 6 

18 

14.85 

50. 7 

19 

15.52 

53. 0 


The heat transfer for test 20 was not presented because only two data points 
were taken which were not enough to determine the slope. These results 
showed that a 14 % change in heat transfer resulted from the presence of the 
thermocouples. 

The percentage of heat transfer attributed to the thermocouples by 
the above methods showed relatively good agreement while the heat transfer 


30 



showed a wide disparity. The heat transfer based on ihe test data was con- 


siderably lower than that based on the analytical solution, and as noted before 
the temperature error went from negative to positive and increased in 
magnitude when the surface coating was added. Based on these results, the 
coating used on the last 10 test runs apparently induced more error (and in 
the opposite direction) than the thermocouples. Further quantitative analysis 
of the thermocouple effect did not appear meaningful because of the coating 
effect and because the heat transfer from the rod cannot be precisely deter- 
mined. The heat input to the specimen was precisely measured but it could 
not be accurately determined what amount of heat was lost from the rod and 
what amount was lost from the heater portion of the test specimen. 

The specimen was coated after the thermocouples were installed and 
therefore the layer of coating actually insulated the surface of the rod and the 
thermocouples. Consequently, the thermocouples were reading the rod 
surface temperature and therefore indicated a temperature higher than the 
surface temperature of the coating. To study this effect, test 19 was used 
because it had the fewest number of thermocouples and yet provided a data 
point at each station. The following numerical differentiation formulae [27] 
were used to calculate the slope of the temperature distribution at each station 
location: 



12AZ + 801 ” 80-1 + 



(IV-1) 

(IV-2) 


i 


[302 + 1001 “ 1800 + 60 - 1 “ 0 - 2] 
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J, 1&2 ~ 60! + 1800 - 100 _i - 30_ 2 ] 


(IV-3) 



-1 


W {2AZ [2502 “ 48e i + 360 O “ 166 -1 + 30 - 2 ] 

2 

» 7^ t“ 30 2 + 160 1 " 360 O + 480 _i - 250 _ 2 ] 


-2 


(IV-4) 

(IV-5) 


Equation (IV-1) was used to find the slope of the curve 3 or more data 
points away from either end of the curve. Equations (IV-3) and (IV-5) were 
used for the 2 points at the beginning of the curve, and equations (IV-2) and 
(IV-4) were used for the 2 points at the end of the curve. The negative 
subscripts denote stations left of zero and the positive subscripts denote 
stations to the right of zero. 

The heat flowing through the rod by conduction at each station was 
then calculated from Fourier's law using the slopes calculated at each station. 
The heat lost by radiation for the rod sections between stations was the 
difference between the heat conducted into the section at the first station and 
the heat conducted out at the next station. This method determined that the 
first section of the rod lost 3. 0 W (10. 2 Btu/hr) and each of the remaining 
sections lost approximately 1. 3 W (4. 3 Btu/hr) . 

The temperature drop across the layer of coating was then calculated 

from Fourier's law assuming that the rod section was iosthermal in the 

longitudinal direction. The carbon black coating was assumed to be 0. 025 cm 

(0. 010 inch) thick and the thermal conductivity value was taken from a table 

in reference 1 which presented data for lamp soot at low pressure. These 
calculations indicated an 18.4*K (33.2*R) temperature drop across the coating 
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on the first section and approximately 7. 8° K ( 14° R) across the coating on the 
remaining sections. 

A corrected temperature distribution was then computed from equation 

(11-18) using an initial temperature equal to the indicated test temperature 

minus 18 . 4° I^\^33 . 2* R) or 333. 0®K (599.37*R) . This curve was plotted in 

Figure C-9 of Appendix C and designated 9 . The test data disagreement was 

c 

reduced from 4. 9 % to 3. 2 % based on the end point temperature. Although 
this did not explain the error in the data, it did indicate a trend in the proper 
direction to explain the positive error for each case where the coating was 
applied. Other variables which could not be accurately determined were the 
actual thickness and the thermal conductivity of the coating. Also, the coating 
was of a porous nature with a velvety surface composed of many small peaks 
and valleys which would radiate from some point beneath the actual surface. 
This compounded the problem of determining the coating thickness. 

Parameter Study 

In addition to the analytical predictions of the test cases, several 

temperature distributions were calculated using equation (11-18) by varying 

the governing dimensionless parameters. The black case was the only one 

considered since for efficient heat transfer the surface would be made as 

black as possible. The configuration factor F^ ^ was plotted versus 

R B 

the length of the rod, in Figure C-10 of Appendix C, for the test specimen 
case to visualize how it affected the temperature distribution equation (IE-18) 
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where 6 






J 


B 


= 1 . 


The above equation was then solved for values of N of 0. 5, 2, and 

c 

10. The dimensionless ratios, Rj and R 2 , which appeared in the term 

F dA - A were varied over a wide range. The limiting cases of no base 
R B 

and an infinite base were also solved for each value of N^. The resulting 
temperature distributions were plotted in Figures C-il, C-12, and C-13 
of Appendix C. 

Some interesting limiting cases of the temperature distribution were 
plotted in Figure C-14. The case of infinite conductivity yields a uniform 
temperature along the rod equal to the base temperature. For the case of 
zero conductivity the governing equation for the temperature of a differential 
element of the rod viewing the base from a distance Z, reduces to 

r 1 i/4 


0(Z) 



(IV-6) 


For a base with an infinite radius (R 2 = °°) , F , A . =0.5 and the tern- 

* » dA - A_ 

R B 

perature along the rod is constant at 9 = (0. 5) 1 / 4 or 0. 8409. Three other 
limiting cases for zero conductivity were also presented in Figure C-14 of 
Appendix C. 

The dimensionless heat transfer or fin efficiency was calculated from 

equation (11-21) for each of the calculated temperature distributions. The 

fin efficiency was plotted versus the conduction parameter for the limiting 

cases of R 2 = Ri (no base) , R 2 = <*> (infinite base, F . A = 0. 5) and 

^R" A B 
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I 


two intermediate values of R 2 for each value of R*. The results are presented 

in Appendix C in Figures C-15, C-16, C-17, C-i8, and C-19. From these 

curves it is seen that a low value of N is required to maintain a high 

c 

efficiency. Also it was observed that the fin efficiency for values of > 0. 05 
was decreased by approximately 50 % from the no base condition to an infinite 
base for all values of considered. 

The temperature distribution curves and efficiency curves are 
applicable to gray surfaces if only singly reflected radiation rays are con- 
sidered. It is only necessary to change the N to N' where N ! = eN . 

c c c c 

CONCLUSIONS 

The experimental data and the numerical solutions showed reasonable 
trend agreement. The experimental data, however, were substantially 
affected by two identified sources of error: 

1. the heat conducted from the rod by the thermocouples 

2. the insulation effect of the carbon coating covering the rod surface 
and the thermocouple junction. 

The above mentioned errors would tend to invalidate the assumption of one- 
dimensional heat transfer. Also the fact that the thermal properties of the 
specimen material were assumed constant would induce some error since e 
and k would vary somewhat with temperature. The emissivity could vary with 
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time from test to test since some surface oxidation would occur. As a result 


of these errors the data did not agree with the analytical predictions as closely 
as was desirable for use in quantitative predictions. The trends did agree, 
however, and the data were within the range of the identified errors in the 
experiment. 

The calculated temperature distributions for a given value of the 
conduction parameter were bounded by the cases of no base, and an infinite 
base; the temperature spread between the two cases increased with increasing 
N . The temperature distribution traversed the range between the limiting 
cases as the values of Rj and R 2 were varied such that the configuration factor 
between the rod and the base varied from 0 to 0. 5. 

The heat transfer from a rod or the rod efficiency may vary a large 
amount with changing base size. The variation may be as much as approxi- 
mately 50% for values of N 0. 05. 

c 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, July 1968 
124-09-05-00-62 
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APPENDIX A. 


CONFIGURATION FACTOR ANALYSIS 


Configuration Factor F s/\ D - a d 

R B 


The configuration factor F of equation (II-18) was determined 

6A r ~ a b 

by the contour integral method [26] . The basic equation for the configuration 
factor for an infinitesimal element on the rod to that portion of the base which 
it views is 


cos cp cos cp _ 

F dA„ . A' - /, 7* ' I, ' (A -‘> 

R B A fi A b 


Applying Stokes theorem to the above area integral results in the following 
contour integral 


R B 


'nf 


y R 7 B 


2tt S 2 

< x B- X R )dZ B- (Z B 

- V^B 

(A-2) 

+ m R 

J 27 T S 2 

* (y B- y R )dX B- (X B 

' X R )dy B 

+ n Ri 

V 2tt S 2 




where 1 , m , n are direction cosines of the normal to the surface dA 
R R R -tv 

and the cosine of the angles, cp and cp , between the normals and the 

B xv 
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vector S, connecting the differential areas has been written as the scalar 

products of n • § ^ and n_ • S_„. The contour integral equation was 
* R RB B BK 

s im plified by orienting the coordinate axes so that the normal to dA R lies 
along the x axes, as shown in Figure A-i. Then by inspection, it was seen 
that 1 R = 1 and m R = n R = 0; also x R = r R , y R = z R = 0, z R = c and dz R = 0. 
The distance S between an infinitesimal area dA R and dA R was 

s = Vu b -v ! + ^b + 4 (A - 3) 


and the contour integral equation (A-2) reduced to 

d y T 


dA 


R 


l B " L 


2tt[ (x. 


R 


B B 

r ) 2 + yL + zi, ] 

r' ^b b j 


(A-4) 
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FIGURE A-i. BASE SURFACE AREA VIEWED BY DIFFERENTIAL 

ELEMENT ON ROD 


The contour that bounded the area A' was composed of two parts as 

B 

shown in Figure A- IB: an arc running from -<p 0 ^ cp ^ cpQ and a straight line 

running from -r^ sincp^ ^ ^ r^ sin^g. On the arc = r^ cos cp and 

y = r sin cp: whereas on the striaght line x = r coscp Qi and dx = 0. 

B B T B ±5 B 

When equation (A-4) was integrated along the arc and then along the straight 
line of Figure A-i. B, it became 


J B 




/ 

<Po 


r cos cp dtp 
B 

r 2 + r 2 + z 2 - 2r r cos <p 


R 


R B 

d y. 


y B 


1b r o 

* -r silled <r B 
B 


(A-5) 


When sin^o and cos<p 0 were written as 


sin<p 0 = 


r 2 - r 2 
B R 

' r_ 


COS<^q — 


JR 

r B 


and r 2 + r 2 + z 2 = C and -2r r = D, equation (A-5) then became 
B ft B Jti B 


^R^B 


Z B r B cos d cp + Z B ^ 

7T j C + D COS cp 7T j 


dy 


B 


<Pa ' r R J B 


~ yt + Z B 


Integration of equation (A-6) and substitution of limits yielded 


(A— 6) 


1 _l 

F . , = — tan 1 

^R " A B n 


1 


r 2 - r 2 
B R 


2Z B Cr B 

7 r D \] C 2 - D 2 


tan 


-l 


(C - D) tan 


Z B r B^° 
7 r D 


( A-7) 
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The term tan was expressed as 


tan 


( ( Pn\_ si 
V 2 ) 1 + 


sin <pn 
cos (f o 


f 


r 2 -r 2 
B R 


r - r_ 
B R 


( A-8) 


and substituted into equation (A-7) . When the appropriate substitutions were 
made for C and D and when <p 0 was written as tan -1 , the expression for the 
configuration factor became 


F,. a, =— tan -1 
^R " A B * 


r 2 - r 2 
B R( 


B 


r 2 + z 2 ) 


B' B 


R 


B 


"rV* 


r B +r R +Z B )2 - 4r B r R 


tan 


27rr 


-l 


l(r B + r R ,J+Z BJl/ r B- r : 


lr B + 'R 1 V (r B +r R +2 B ,! ‘ 4r B r : 


■2 ^2 
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tan 
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R 


r 2 -r 2 
B R 


R 


( A-9) 


When the above equation was nondimensionalized with respect to L, and R 1( 

R 2 , and Z were defined as Rj = r /L, R 2 = r /L, and Z = z /L, the 

R B B 

nondimensional configuration factor became 


dA - A' =— tan -1 — R * 

R B 7 T f 


2 ~ 

z 


z(R 2 + R 2 + z 2 ) 


tan 


-l 


) ttR* Y(R§+ Rf + z ^) 2 - 4Ri Rf 
[ (R 2 + Rj) 2 + Z 2 ] Yr \ - Rj 2 


(R 2 + Rj) V(Rl + R? + z 2 ) 2 - 4R 2 R 2 


2ttRi 


tan" 


R z -qZ 

2 - Rl 

Ri 
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and from symmetry 


F 





(A-ii) 


Configuration Factors F, Ar . sAg and F s Ag . 8Ar 


It was readily seen from the cylindrical symmetry of the specimen 
that the flux density at any arbitrary radius p on the base was constant for any 
arbitrary angle cp. The same was true for any cp position on the rod at a 
given z. 

Now when a differential slice of the rod and a differential ring on the 
base were considered, it was seen that the radiosity distribution was uniform 
with the coordinate angle cp . However, the radiosity distribution varied along 
the length of the rod and also along any radius on the base. 

Also from symmetry, the configuration factor from a differential 
element on the rod to the portion of a differential ring viewed on the base was 
equal to the configuration factor from the total differential ring around the 
rod to the total differential ring on the base, or 



- 6A[ 


B 




( A-12) 


then 


c^R-^B 


dA 


R 


/ £ <* 


6A k 


B’ 


«, R> S) dA R dA B , 


(A-13) 
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or 

F 6A - SA = f f ^B ,V R*^ ‘ (A-14) 

R B &A' 

B 

The differential areas were expressed as dA = r dzd^" and dA = pdpdcp, 

XV R B 

as shown in Figure A -2. 

y 



FIGURE A-2. COORDINATE ANGLE THROUGH WHICH EITHER 

dA^ CAN VIEW <5 A OR dA CAN VIEW 6A n 
K B B R 

The total coordinate angle through which either differential element dA 

R 

and dAg could view the differential strips 6A-g or <5A R was 2cos ^r-^/p. 
42 


Therefore, the configuration factor from the differential strip 6A to the 

Jet 


differential strip 6A could be written as 2 times the integral from 0 to 
B 

-i r R 

cos 1 or, 

P 

cos 1 _ 

F 6A„-aA - 2 / P <“«• 

R B 0 


B 


(A-15) 


From the reciprocity relationship, 

6A R „ 


6A b ~ 6A r 6A b 6A r- 6A b 


27rr dz 

— F 

27rpdp 6A - 6A 

K B 


, r R dz 

6A b “ 6A r pdp 


-tl* 


cos 

2 f 

o 




(A-16) 


The vector S was written as 


A A A A 

S RB =(X - r R } i + yj " Zk 


( A-i7) 


and 


a a a a 

S BR = ( r R ~ x) 1 - yj + z k 


(A-18) 


and the magnitude of vector S was 


s = V<x - r R ) 2 + y 2 + z 2 . 


(A-19) 


The normals to the rod and base areas were 


a a 


A A A 


n R = i + ( °) j + ( °) k = i 


( A-20) 


43 



( A-21) 


A A AAA 

n B = ( 0) i + ( 0) j + k = k 


The cosines of the angle between the normals and the vector S were found by 

A 

taking the vector dot product of the normal and the vector S and dividing by 
the product of the magnitude of £> and the magnitude of the unit vector, which 
yielded 


x - r 


COS <p R = 


R 


V (x - v 


2 + y 2 + z 2 


(A-22) 


and 


cos <p B = 


V(x - r R ) 2 + y 2 + z 2 


( A-23) 


These were converted to cylindrical coordinates by substitution of x = pcoscp 
and y = pshkp which yielded 

pcoscp - r 

Jt\ 


cos <p R = 


y p2 + r R + z2 ■ 2r R pCOS< ? 


( A-24) 


cos cp B = 


yV + r^ + z 2 - 2r R pcos^T 


( A-25) 


S 2 = p 2 + r 2 + z 2 - 2r pcoscp . 
R rt 


(A-26) 


It was then written that 


( pcoscp - r ) ( z) 


7r(p 2 + r^ + z 2 - 2r R pcos <p) 2_ 


R' 


(A-27) 


and substitution into equations (A-15) and (A-16) yielded 
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,-i_R 


6A r - 6A b 


- COS ” 1 (pCOSrt? - r p) 

Pdp r P # 5 7 

7 r J (p 2 + r 2 + z 2 - 2r_ pcoscp)' 

Q XV XV 


„ i ±t 

2r zdz cos - 

R r P 


6A b - 6A r 


( pcoscp - r R ) 

( p 2 + r 2 + z 2 - 2r pcosm) ‘ 
xv It 


When the integral in equations (A-28) and (A-29) was considered, 


cos 1 (pcoscp - r^) 

f (p 2 + r 2 + z 2 - 2r pcoscp ) T ^ 

0 K K 


The integral (A-30) became 




cos 1 (pcoscp - r ) < 

f p 7 £-a cp= f 

J ( a + b cosep) r J 


p Ad<p 

(a + b cos<p) 


Bd<p 

(a + b cos<p) 


when a and b were defined as a = p 2 + r 2 + z 2 and b = -2r p and the 

XV XV 

expression under the integral sign in equation (A-30) was separated, using 
partial fractions, 


p G° s <p - r _ A B 

(a + b cos cp ) 2 (a + b cos cp) 2 (a+beos cp) 5 


where B = — and A = 
b 


pa 

r + — — 
R b 


45 



From the integral tables [28] the above integral was found to be 


/ 


R 

cos -1 ( pcoscp - T )d(p 

p K T- = M 


(a + b cos cp) ‘ 


( b 2 - a 2 ) ( a + b cosep) 


b sin cp 


~ b^r? 


tan -1 


( a - b) tan 


(f) 






2B A 
+ — tan 


( a - b) tan 


(i) 


cos 


L R 


1 


Va 2 - b 2 


( A-33) 


When the integration limits and the following relationships 


cosep 


R 




, sin<p =■ 


R 


tan 


(D-i? 


sin^? 




R 


co scp p + r 


R 


were substituted into equation (A-33) the expression became 


/ 


cos 


-1 


R 


{peosep - r R )d<^ 


A b 


V ? 7 


R 


(a + b cos cp) ‘ 


p(b 2 - a 2 )^a+ b — : 


2aA 


(b 2 - a 2 ) Va 2 - b 2 


.tan 


-i 


(a - b) yp “ r 


4 


( P + r R ) V a 2 - b 2 


2B , _i 
+ — ■ ■ tan 


V7T 


(a - b) Yp 2 - r: 
(p + r R )V a 2 - b 2 


(A-34) 


and this was further simplified by substituting for A and B and combining 
terms to yield 
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/ 


cos 


.ii 

p 


( pcoscp - r R )d<p yp 2 - r: 
T~ ="T? 


(a + b cos<p) ‘ 


ar - b* 


2(ar + pb) 

5 r— tan" 1 

(a 2 - b 2 ) 3/2 


(a-b) (p- r R ) 

( a + b) (p + r R ) j 


(A-35) 


When the above expression was used and appropriate substitutions were made 
for a and b, equation (A-28) for the configuration factor became 


m r - “b 


_ 2zpdp 
7T 


'/p 2 - r 2 


R 
2 \ 2 


[(p 2 + r^+ z 2 ) 2 - 4p 2 r R ] 


2r R (p 2 - - z 2 ) 


[(p z + r* + z 2 ) 2 - 4r*p £ ] 


„2 J. 


,3/2 


R 


R 


-1 


tan 


Similarly, equation (A-29) became 


(p 2 + r R + z 2 + 2r R p) (p - r R ) 
(p 2 + r^ + z 2 - 2r R p) (p + r R ) 


(A-36) 


2r zdz 
R 


*b-“r 


I(p r + + z 2 ) - 4 P 2 r 2 n ] 


R J 


2r R (p _r R “ Z ' 




tan 


-l 


(p 2 + r^ + z 2 + 2r R p) (p - r R ) 
(p 2 + r^+ z 2 - 2r R p)(p+ r R ) 


( A-37) 


Equations (A-36) and (A-37) were expressed in the following form: 


47 



when the dimensionless ratios Z = z/L, Rj = r^/L, R = p/L were introduced. 


48 


APPENDIX B 


EXPERIMENTAL EQUIPMENT 


The experimental equipment consisted of the high altitude simulation 
system, the data acquisition system, and the emissivity measuring device. 

This equipment is government property maintained and operated by the 
Propulsion and Vehicle Engineering Laboratory of the George C. Marshall 
Space Flight Center. 

High Altitude Simulation System 

The high altitude simulation system or vacuum system consisted of the 
following primary components: vacuum chamber; diffusion pump, ejector pump, 
backing pump, and roughing pump; valves; pressure gauges; and liquid nitrogen 
storage trailer. 

Vacuum Chamber. The vacuum chamber was a stainless steel cylinder 
approximately 1. 83 m (6 feet) long and 1.22 m (4 feet) in diameter with a 
dished head at each end. One of the dished heads was hinged for a door and the 
other attached to the pumping system. The chamber environment was controlled 
by a cooling jacket which permitted a cooling liquid such as liquid nitrogen to be 
circulated through the cooling jacket and the chamber door. The door was a 
double wall construction connected to the jacket so that the cooling liquid 
circulated through both the jacket and the door. The interior surfaces of the 
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chamber that were viewed by the test specimen were painted black to minimize 
radiant energy reflections. 

The vacuum chamber wall was designed with 3 ports that could be used 
for viewing or for instrumentation or for power feedthrough. For this test 
the ports were used for plugs to feedthrough 64 thermocouples and 2 power 
supplies. 

Pumping System. The pumping system consisted of 4 pumps which could 
produce a vacuum chamber minimum pressure of 1 x 10“ 7 mm Hg. The 
chamber vacuum was initiated by a roughing pump of large capacity. After a 
rough vacuum was obtained (150 mm Hg) , the low vacuum was achieved and 
maintained by a fractionating oil diffusion pump, an oil ejector pump, and a 
rotary backing pump operating in series. The pumping system was designed 
and manufactured by the Consolidated Electro-dynamics Corporation for NASA/ 
MSFC. The following pumps were used in the system: 

1. Type MCF-15000 fractionating oil diffusion pump with 0. 81-m 
(32-inch) diameter casing. 

Operating Range: 5 x 10 -7 to 5 x 10 -4 mm Hg 

Maximum Speed: 899. 9 m 3 /man (31, 779 cfm) at 3 x 10 -4 mm Hg 

2. Kinney Type KB-150 oil ejector pump. 

Operating Range: 6 x 10 -4 to 4 x 10 -1 mm Hg 

Maximum Speed: 11. 53 mVmrn (407 cfm) at 3 x 10 -2 mm Hg 
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3. Kinney Type KDH-130 rotary oil sealed backing pump. 

Free Air Displacement: 3.71 m 3 /min (131 cfm) at 535 rpm 
Ultimate Pressure: 1 x 10 -2 mm Hg 

4. Kinney Type KD-310 rotary oil sealed roughing pump. 

Free Air Displacement: 8. 81 m 3 / min (311 cfm) at 360 rpm 
Ultimate Pressure: 1 x 10 -2 mm Hg 
Valves and Gauges. The high altitude simulation system design 
incorporated the necessary valves and pressure gauges to control the system. 
All system valves were hand operated except the 0. 81-m (32-inch) valve 
between the vacuum chamber and the diffusion pump which was motor driven. 
There were 3 types of pressure gauges used in the system: a cold cathode 
ionization and a Pirani type were used to monitor the chamber vacuum; a 
Bourdon-tube gauge was used to monitor the pressure in the liquid nitrogen 
storage trailer. Both the high and low vacuum gauges were products of the 
Consolidated Vacuum Corporation, Rochester Division, Rochester, New York, 
and were of the following specifications: 

1. Cold cathode ionization type 

Discharge Vacuum Gauge 
TYPE GPH-100A 

Three Range Selections from 25 x 10 -3 to 1 x 10 -8 millimeters 
of mercury 
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2. Pirani Vacuum Gauge 


TYPE 2203-04 

Two range selection 2000-50 microns 
The Pirani gauge was used to monitor the vacuum while the system was 
operating at a low vacuum and the ionization gauge was used to monitor the 
vacuum when the system was operating at a high vacuum. 

Liquid Nitroge n Storage Trailer. Liquid nitrogen was required for the 
vacuum pumping system cold trap and also for the vacuum chamber cooling 
jacket. The liquid nitrogen was supplied from an 8165. 0-kg (9-ton) capacity 
liquid nitrogen storage trailer. The storage trailer was a self-pressurizing unit, 
and was maintained at a pressure of 1. 72 x 10 5 - 2 . 07 x 10 5 N/m 2 (25-30 psig) 
for the test runs. Flow was controlled to the vacuum system by a hand valve and 
the storage trailer pressure was maintained by a hand vent valve and an auto- 
matic (spring-loaded) relief valve. The storage trailer was located outside of 
the building which houses the vacuum system for safety reasons, and the liquid 
nitrogen was supplied by a 2. 54-cm (1-inch) diameter insulated pipe. 

Data Acquisitio n System. The data acquisition system consisted of 64 
thermocouples, a digital recording system, and a tape controlled typewriter. 

The thermocouples were standard 30-gauge copper and constantan 
thermocouple wiring. The thermocouples were spot welded to the specimen 
and brought out of the vacuum chamber through a standard vacuum type feed- 
through electrical connector. A thermocouple cold junction (ice bath) was 
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located adjacent to the vacuum chamber and the thermocouple outputs were 
carried by way of a floor conduit to the recording system approximately 
15.24 m (50 feet) away. 

The digital recording system was a model ADRT 2100 manufactured 
by the Minneapolis-Honeywell Regulator Company. The system could record 
up to 100 channels of data. The system could successively sample and record 
data inputs at the rate of 100 points per minute through the use of stepping 
switches. A sampled input was amplified and then balanced against a precision 
voltage. The balance was obtained by successive approximation performed by 
relays pulling in successively decreasing increments of voltage. If an increment 
pulled in was too large, it was dropped and the next smaller increment pulled in; 
if the latter was too small, another smaller increment was pulled in, etc. , , 
until the balance was achieved. A digital output was then punched on paper tape 
along with the necessary control signals to control a Friden Flexowriter. The 
tape was fed through the Flexowriter control system at each recording interval 
and the data outputs were typed out in millivolt units. The data acquisition 
system could scan the selected number of inputs and punch the output values on 
paper tape once every 1/2, 1, 2, 5, 10, 20, 50, or 100 minutes, or continuously. 

For this series of tests, 64 thermocouples were monitored and recorded 
at intervals of 5 minutes from test beginning to completion. All thermocouple 
millivolt outputs were available for individual call up at any time. In the call 
up mode, a digital reading appeared in lights by making the proper manual 
selection to read the particular measurement of interest. 
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The data acquisition system accuracy was dependent upon the range 
selected. For this test series a range of 0 to 10 millivolts was selected and 
the system accuracy in converting the thermocouple output was ±0. 01 millivolts 
[approximately 0. 22° K (0. 4° F) ] . 

Emissivity Measuring Device 

The emissivity of the base surface specimen was measured by a Gier 
Dunkle-Emissivity Inspection System, Model EM 520, Geir Dunkle Instruments, 
Inc. , Santa Monica, California. The emissivity inspection system measured 
the emittance of an opaque surface at room temperature. The system was 
composed of a heated radiometer connecting a flexible conduit to a carrying 
case containing the power supply. When a room temperature specimen surface 
covered the radiometer opening, the detector provided an output that was 
proportional to the infrared emittance of the specimen. The output was 
indicated on the self balancing potentiometer and this output was then converted 
to an emissivity valve. The emissivity measurement of the base surface 
specimen was made by personnel of Materials Division of Propulsion and 
Vehicle Engineering Laboratory and the emissivity was found to be 0. 13. 
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APPENDIX C 


ILLUSTRATIONS 
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FIGURE C-3 


SCHEMATIC OF TEST SETUP 













TEMPERATURE RATIO, 0 = T/T 



FIGURE C-4. TESTS i-4 AND CALCULATED TEMPERATURE 
DISTRIBUTIONS FOR ROD WITH BASE (Uncoated) 


TEMPERATURE RATIO, 0 = T/T, 
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o 0.5 1.0 


LENGTH RATIO, Z = z/L 

FIGURE C-6. TESTS 8-13 AND CALCULATED TEMPERATURE 
DISTRIBUTIONS FOR ROD WITH NO BASE ( Coated and Uncoated) 


LENGTH RATIO, Z = z/L 


FIGURE C-7. TESTS 14-17 AND CALCULATED TEMPERATURE 
DISTRIBUTIONS FOR ROD WITH BASE ( Coated) 




310 


OTEST NO. 12 
□ TEST NO. 18 
A TEST NO. 19 
OTEST NO. 20 


0 



LENGTH RATIO, Z = z/L 


FIGURE C-8. TEMPERATURE DISTRIBUTIONS FOR THERMOCOUPLE 
CORRECTION TEST RUNS, T(*K) (Coated) 





LENGTH RATIO, Z = z/L 


FIGURE C-9. TEMPERATURE DISTRIBUTIONS FOR THERMOCOUPLE 
CORRECTION TEST RUNS, 0 = T/T 0 (Coated) 



CONFIGURATION FACTOR, F dA 


I 



o 10 20 30 40 50 

DISTANCE, z ( CENTIMETERS ) 


05 

U1 


FIGURE C-10. CONFIGURATION FACTOR OF A DIFFERENTIAL ELEMENT 
OF ROD VIEWING RASE FROM DISTANCE z 
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FIGURE C-li. TEMPERATURE DISTRIBUTIONS FOR FIXED 
N AND VARIOUS Rj AND R 2 RATIOS ( Black Surfaces) 



TEMPERATURE RATIO, 0= T/T ( 



FIGURE C-12. TEMPERATURE DISTRIBUTION FOR FIXED 
N c AND VARIOUS R t AND R 2 RATIOS (Black Surfaces) 
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LENGTH RATIO, Z = i L 

FIGURE C-13. TEMPERATURE DISTRIBUTIONS FOR FIXED 
N c AND VARIOUS R t AND R 2 RATIOS (Black Surfaces) 
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TEMPERATURE RATIO 
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LENGTH RATIO, Z = z/L 



FIGURE C-14. SOME LIMITING CASES OF TEMPERATURE 
DISTRIBUTION FOR BLACK SURFACES 
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FIN EFFICIENCY 



FIGURE C-16. FIN EFFICIENCY VARIATION WITH CONDUCTION PARAMETER FOR FIXED 
VALUES OF Rj AND VARIOUS R 2 RATIOS, R t = 0. 05 (Black Surfaces) 














TABLE C-i. CALCULATED TEMPERATURE DISTRIBUTIONS 


Dimensionless 
Length, Z 

Conduction Parameter, ‘N = 

c 

2.0 

- 0. i and R 2 = 0 

Rj = 0. 1 and R 2 = 2. 0 

Ri = 0. 1 and R 2 = 00 

Dimensionless 
Temperature, q 

Dimensionless 
Temperature, 0 

Dimensionless 
Temperature, 0 

0.0 

1.00000 

1.00000 

1.00000 

0. 1 

0. 92712 

0. 95783 

0. 96708 

0. 2 

0. 86916 

0. 92324 

0. 94171 

0.3 

0.82271 

0. 89454 

0. 92212 

0.4 

0. 78549 

0. 87064 

0.90701 

0. 5 

0. 75591 

0. 85083 

0.89546 

0.6 

0. 73289 

0. 83468 

0. 88678 

0.7 

0. 71566 

0. 82201 

0. 88048 

© 

00 

0. 70369 

0. 81274 

0.87621 

0. 9 

0. 69664 

0. 80699 

0. 87374 

1.0 

0. 69431 

0. 80496 

0. 87293 
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TABLE C-2. TEST SUMMARY 


Test 

No. 

Date 

Rod 

Power (W) 

Base 

Power (W) 

Vacuum 
(mm of Hg) 

Test 

Duration 

(Hours) 

1 

6/16/67 

40.4 

0 

1.5 x 10~ 5 

6.5 

2 

7/11/67 

40.4 

0 

4.0 x 10~ 7 

8. 0 

3 

7/17/67 

40.4 

30. 6 

1.2 x 10 -6 

8. 0 

4 

7/19/67 

40.4 

42.2 

2.0 x 10" 6 

8. 0 

5 

7/20/67 

18. 6 

0 

2. 5 x 10~ 6 

8. 0 

6 

7/21/67 

18. 6 

22. 5 

3. 0 x 10“ 6 

8. 0 

7 

7/24/67 

18. 6 

10. 0 

1.2 x 10- 6 

8. 0 

8 

7/25/67 

18. 6 

No Base 

1. 5 x 10 -6 

9. 0 

9 

7/25/67 

29. 0 

No Base 

8. 0 x 10- 7 

12. 6 

10 

7/26/67 

40.4 

No Base 

8. 0 x 10 -7 

8.8 

11 

8/7/67 a 

18. 1 

No Base 

7.0 x 10 -7 

13. 7 

12 

8/7/67 a 

41. 6 

No Base 

5.0 x 10" 7 

10.5 

13 

8/8/67 a 

75. 3 

No Base 

3.0 x 10 -7 

11. 5 

14 

8/9/67 a 

41. 8 

10. 0 

5. 0 x 10“ 7 

17. 3 

15 

8/10/67 a 

41. 6 

62. 5 

8. 0 x 10" 7 

11. 0 

16 

8/10/67 a 

41. 6 

90. 0 

8. 0 x 10 -7 

11. 8 

17 

8/ll/67 a 

40.4 

22. 5 

3. 0 x 10' 7 

10. 8 

18 

8/15/67 a 

41. 6 

No Base 

1.2 x 10“ 6 

11.2 

19 

8/17/67 a 

41. 6 

No Base 

1. 5 x 10~ 6 

10. 9 

20 

8/18/67 a 

41. 6 

No Base 

9.0 x 10- 7 

11. 5 


a. Denotes the specimen was carbon coated. 



TABLE C-3. TEST DATA 


Test 

Number 

Temp, of 
Base (T ~ °K) 

.D 

Temp, of Rod 
At Base (T 0 ~ °K) 

N 

c 

1 

393.0 

508. 6 a 

1. 70488 

2 

365.2 

457.4 

1. 24091 

3 

477.2 

497. 6 

1. 59716 

4 

512. 9 

503.4 

1. 65317 

5 

305. 2 

353. 5 

0. 57264 I 

6 

427. 8 

380. 5 a 

0. 71391 

7 

407. 3 

380. 5 

0. 71391 

8 

— 

343. 5 

0. 52541 

9 

— 

389. 9 

0. 76818 

10 

— 

445. 4 

1. 14507 

11 

— 

280. 6 

0. 28617 

12 

— ! 

345. 8 

0. 53601 

13 

i 

423. 6 

0. 98544 

14 

344. 6 

368.7 

0. 64961 

15 

449. 9 

394.4 

0. 79549 

16 

490. 2 

403. 2 

0. 84957 

17 

367. 2 

377. 2 

0. 69586 

18 

— 

348. 8 

0. 54987 

19 

— 

351.4 

0. 56258 

20 

— 

352. 1 

0. 56579 


a. Steady state conditions were marginal. 
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TABLE C-4. TEMPERATURE ERROR AT THE INSULATED 

END OF THE ROD 


Test 

Number 

Dimensionless 
Temperature, 0 u, 
of Rod ( Test) 

1 

0. 9113 

2 

0. 9159 

3 

0. 9003 

4 

0.8988 

5 

0. 9525 

6 

0. 9381 

7 

0. 9430 

8 

0. 9560 

9 

0. 9411 

10 

0. 9210 

11 

0. 9344 

12 

0.8883 

13 

0. 8322 

14 

0.8717 

15 

0. 8568 

16 

0. 8512 

17 

0.8668 

18 

0. 8829 

19 

0. 8822 

20 

0. 8840 




Dimensionless 
Temperature, 6 u, 
of Rod ( Calculated) 


0. 9169 
0. 9354 
0. 9210 
0. 9189 
0. 9667 
0. 9596 
0. 9596 
0. 9719 
0. 9597 
0. 9431 
0. 8960 
0. 8473 
0. 7865 
0. 8363 
0. 8181 
0. 8160 
0. 8323 
0. 8433 
0. 8408 
0. 8396 
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